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ABSTRACT 

The three dimensional equations of motion for a cable connected 
space station — counterweight system are developed using a Lagrangian 
formulation. The system model employed allows for cable and end body 
damping and restoring effects. The equations are then linearized about 
the equilibrium motion and nondimensionalized. 

To first degree, the out-of-plane equations uncouple from the in- 
plane equations. Therefore, the characteristic polynomials for the 
in-plane and out-of-plane equations are developed and treated separate- 
ly. From the general in-plane characteristic equation, necessary con- 
ditions for stability are obtained. The Routh-Hurwi tz necessary and 
sufficient conditions for stability are derived for the general out-of- 
plane characteristic equation. Special cases of the in-plane and out- 
of-plane equations (such as identical end masses, and when the cable is 
attached to the centers of mass of the two end bodies) are then examined 
for stability criteria. 

Time constants for the least damped mode are obtained for a range 
of system parameters by numerical examination of the roots of the in- 
plane and out-of-plane characteristic polynomials. For the in-plane 
case, a comparison with results previously obtained in a two dimensional 
treatment (but with a different damping scheme) is made. 

The effect of first order gravi ty- gradient torques on the steady- 
state motion is shown to be small. Resonance due to gravity-gradient 
forcing terms is examined and is seen to occur for certain choices of 
system parameters. 
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I. INTRODUCTION 

Artificial gravity in a space station system may be created in 
two different ways: the station in a rim-like configuration may be 

rotated about its axis of symmetry, or the station connected to a 
counterweight by a taut cable can be rotated about the system center 
of mass. The second system may have certain weight and power system 
advantages over the first; to change the spin-rate of the system it 
is necessary to adjust the effective equilibrium length of the cable, 
whereas in the rim configuration an active power source is required. 

One of the earlier treatment studying the dynamic behavior of 

cable-connected two-body systems was given by Paul who considered 

the planar motion and stability of a gravity-gradient stabilized, 

extensible dumbbell satellite system where the cable mass effects 

were neglected. Paul developed stability criteria and showed that, 

if the internal friction resulted from "material damping" within the 

elastic cable, there would be relatively little damping of a viscous 

nature, but that a nonlinear time-independent type of hysteretic 

2 

damping could be significant. Bainum et al included the effects of 
distributed (unsymmetri cal ) end masses for the case of a gravitation- 
ally stabilized tethered-connected interferometer system and con- 
cluded that a combination of tether system damping and rotational 
damping of the motion of the end masses about their own mass centers 
must be employed; the use of one damping scheme without the other 
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will not provide adequate damping of each normal mode for the rela- 
tively long (e.g. 3 n. mi.) connecting lengths required. The first 

three-dimensional stability analysis of a tether-connected gravita- 

3 

tionally stabilized system was presented by Robe. His system con- 
sisted of two identical but unsymmetrical distributed end masses con- 
nected by a massless, extensible tether, resulting in nine degrees- 
of- freedom. It was shown that there is a decoupling of the small - 
amplitude motions within the orbital plane from those outside the 
plane; therefore, additional "out-of-plane 11 stability criteria, would, 
in general, have to be satisfied in comparison with the previous two 
dimensional criteria. A paper by Beletskii and Novikova considered 
domains of possible three dimensional motion for a gravitationally 
stabilized point-mass system connected by a flexible, massless tether 
for the cases of both a taut and a slack tether. 

In the area of rotating connected systems an earlier paper by 
5 

Chobotov included the effects of cable mass and elasticity with 
point mass end masses and two-dimensional motion. It was found that 
the gravity-gradient effects upon the small amplitude vibration 
stability of the rotating system are very small and that the stabil- 
ity criteria are functions of the cable natural frequencies, the 

angular velocities of the station and orbital motion, and viscous 

6 

damping parameters. Subsequently, Stabekis and Bainum examined the 

motion and stability of a rotating space station-massless cable- 

( 

counterweight configuration where the motion was restricted to the 
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orbita/ plane. Although the system remained stable in the absence of 
rotational damping (of end body motions), this damping in addition to 
cable damping is necessary to achieve reasonable time constants for 
the nominal parameters considered. A paper by Nixon deals with 
determining the dynamic equilibrium states in three dimensions for a 
completely' undamped system with an arbitrary number of cables. He 
showed that his linear model accurately compares to the nonlinear 
model when the cable tension has some initial value and when angles 
do not deviate from their equilibrium values by more than one or two 
degrees. Although Nixon analytically determined the states of equi- 
librium in three dimensions, he did not perform an analytical stabil- 

8 

ity analysis about these motions. Anderson, whose system had dis- 
tributed end masses with lateral oscillations for three dimensional 
motion, used an energy approach to analyze the motion of the system 
under the influence of disturbance torques. He found that the basic 
attitude response of the space station is that of an undamped second- 
order system and that coupled to this basic response are rigid body 
characteristics and cable lateral mode effects. 

Of interest in this investigation is an examination of the 
three-dimensional motion of the rotating cable-connected system for 
the general case where the end bodies have a distributed mass (finite, 
unequal moments of inertia) and the possibility of energy dissipation 
in both the cable system and end bodies is included. To date* this 
treatment has not appeared in the open literature and would represent 


4 . 


an extension to the problem considered in Ref. 6. Even if a decou- 
pling of the small -amplitude in-plane motions from the out-of-plane 
motions would result (as in Ref. 3), the additional stability cri- 
teria emanating from these out-of-plane motions would have to be 
carefully considered prior to the design of such a system as a means 
of creating artificial gravity. The optimum design parameters of 
such a control system may vary considerably from those inferred from 
a two dimensional analysis. 
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/ 

II. DESCRIPTION of mathematical model 

It is assumed that the system center of mass follows a circular 
orbit, that the cable is extensible but massless and that the system 
at equilibrium has a nominal spin rate in the orbit plane about an 

axis passing through its center of mass. 

Five different coordinate systems describe the motion. The 
fixed inertial reference is located at the center of mass of the 
Earth, whereas, the A coordinate system is located at the center of 
mass of the system model with the A 1 axis along the local vertical, 
the A2 axis in the direction of the velocity of the orbit and the A3 
axis normal to the orbit plane. The third coordinate system is the 
B system fixed in the space station (body 1 as shown in Fig. 1) at 
its center of mass. The axes of the B system are assumed to be the 
principal axes of -body 1 with the cable attached at a point on the 
B] axis. A one-two-three sequence of rotations, respectively, is^ 
assumed to orient the B system with respect to the A system. The C 
system fixed in a body 2 (the counterweight) at its center of mass 
is defined the same way as the B system. Lastly, there is the D 
coordinate system which is located at the center of mass of the 
model and is defined by two rotations with respect to the A system: 
an angle e-j in the orbit plane and then an angle e 2 out of the plane 
By these rotations, the D-j axis is parallel to the cable line. 
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The transformations from the A to the B systems, from the A to 
the C systems, and from the A to the D systems are given in equa- 

tions (1) and (2) as: 



f 


V 


COS 6 3 

sing 3 

“1 

0 


COSB2 

0 

-sir\s 2 

1 

0 

0 


Ai 

b 2 

= 

-sin&3 

COS$3 

0 


0 

1 

0 

0 

cos$i 

sinBi 


a 2 

b 3 


0 

0 

1 


sin $2 

0 

COSB2 

0 

-sinBi 

cosBi 


_ A 3 _ 


Cl 


— 

COS T3 

siny 3 

0 


C0Sy 2 

0 

. _ 

-sin Y 2 

1 

0 

0 


Ai 

C 2 

= 

-Siny 3 

cosy 3 

0 


0 

1 

0 

0 

. COS Yi 

sinyi 


a 2 

C 3 

L — 


0 

0 

1 

J 


siny 2 

0 

COS Y 2_ 

0 

-sin Y1 

COSYi _ 



Di 


COS02 

0 

sin 02 


COS01 

sin 0 i 

0 

Ai 


') 


0 2 

s 

0 

1 

0 


-sin© 1 

COS01 

0 

a 2 

■ 



Ds 

_ — 


-sin 02 

0 

COS0 2 


0 

0 

1 

A 3 





(la) 


Ob) 


( 2 ) 
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III. DEVELOPMENT OF THE EQUATIONS OF MOTION 


The following definitions of vectors are used throughout this 
section (see Fig. 2).- 


i/o 


Position vector from the center of the Earth to 
thp npnter of mass of body i (i = 1 »2) 


7 Vector from the c.m. of the system to the c.m. of 

i /A 

body i (i = 1 »2) 

l ■ Vector from the attachment point of body 2 to the 

.attachment point of body 1 



r l/P 


r 

2/P 

/ 1— 
A/F 

CO 

B/A 

W B,C/F 


Vector from the c.m. of the Earth to the c.m. of 
the model 

Vector from the attachment point of body 1 to the 
c.m. of body 1 

Vector from the attachment point of body 2 to the 
c.m. of body 2 

Angular velocity of the A coordinate system with 
respect to the fixed inertial reference (F) 
Angular velocity of the B system with respect to 
the A system 

Angular velocity of the B or C systems with res- 
pect to the F system 



Angular velocity of the D system with respect to 
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the F system 

(*) The time derivative of a vector with respect to the 

fixed, inertial reference 

( g q q The time derivative of a vector with respect to 

the noninertial A, B, C, or D systems, respectively 

III. A. Energy Expressions and the Rayleigh Dissipation Function 
To use Lagrange's equations, it is necessary to express the 
total kinetic and potential energy of the system in terms of the 
variables which describe the system's motion. This can be performed 
by means of vector equations whose components are functions of the 
variables. 

The equation for the translational kinetic energy is, 

T T " ^rVl/o + j m 2 r 2/o (3) 

*2 

where r^ 0 = r\j/ 0 : r^ 0 for i = 1,2. From the fact that the A sys- 
tem has a noninertial rotation, equations (4) and (5) may be obtained. 

^l/o^F = ^ r l/A^F +a, A/F x r A/o 

w • 

f r 2/oV = + "’A/F x r A/o ( 5 ) 

Equation (6), arising from the geometry of the system (Fig. 2), is 

r 1/P ‘ r l/A + r 2/A ' r 2/p = 0 ( 6 ) 

From the fact that point A is the system c.m. , 



The model used has two end masses. This indicates that for gen- 
eral unconstrained motion twelve variables are needed. If in addi- 
tion, the variables and l are added for the cable line, there 

will be fifteen variables in all. Equations (6) and (7) can be used 
to express r ]/A and r 2/A in terms of the other vectors. Since the x, 
y, 2 coordinates of each end mass are then eliminated, equations (6) 
and (7) can be considered as six equations of constraint. Thus nine 
independent variables describe the motion of the model and are 
selected as follows: e-j and e 2 for the orientation of the cable with 

respect to the A system: & for the variable cable length; 3-j * $2 S ^3 

for the orientation of body 1 and y*{ > ^ 2 * ^3 5 the orientation of 
body 2. 

Equations (6) and (7) can be combined in order to express the 
vectors r^ A and r 2 / A as: 

F 1/A = [T + "1/P ' F 2/P ] (8) 

r 2/A= - I> + n/P - r 2/p] (9) 

The derivative of equations (8) and (9) are taken, noting that 

u = (&)p + ^ / p ^ F ~ ^ r 2/P^F (1^) 

Then after substituting the results of this differentiation into 
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equations (4) and (5), the following is obtained 

% m 2 u ^ — 

(ri/Jr = — + Wfl/r x 


1/o'F 


m]+m2 


A/F * A/o 


OD 


' . m]U - 

(i Wf = ' + “A/F x r A/o 


( 12 ) 


It should be noted also that 

« • • » 

“ = - ^2/a'f = (n/o)F - ^2/oV 


( 13 ) 


Equation (10) may be further expanded using the general equation 
relating the derivative of a vector in an inertial system to its 
derivative in a noninertial system, 


(£)f = Wq + 


5 d/f 


x l 


04 ) 


( r !/p)p - ( r l/p)B + “B/F x r l/P 


( 15 ) 


( r 2/ p) F ^2/P^C + “C/F x r 2/P 

Since the cable can stretch, {£) = £d^. Furthermore, 


06 ) 


(r, /n L “ (r 0 / J r = 0 since r. . and r 0/ are constant vectors in the 
B and C systems, respectively. 

Equations (11) and (12) when substituted into equation (3) 

give 

_ 1 


] 1 __ ___ 

T T = 2 u (u-u) + j (m 1 +m 2 ) |ioa/f x r/\/ 0 1 2 


( 17 ) 


where the second term is constant for a circular orbit and constant 
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orbital angular velocity. 

The rotational kinetic energy has the form 


T R = + 


08 ) 


since the axes of the B and C systems are the principal axes of 
bodies 1 and 2, respectively. 

The potential energy is, allowing for restoring forces in the 
cable and end bodies, 

GM (m-i+mo) 2 

v = . e 1 + k 

r A/o 



(19) 


where G is the gravitational constant, M e is the mass of the Earth and 
the bars indicate variational variables defined later. The interac- 
tion of on-board magnetics with the ambient magnetic field could pro- 
vide rotational restoring torques derivable from the type of potential 
assumed above and discussed in Ref. 2. 

It is further assumed in connection with (19) that the restoring 
force in the cable is proportional to the linear extension from the 
unstretched length— i.e. the cable obeys Hooke's Law. Damping forces 
are assumed to act through a dissipation function of the form: 


F-ll + k Ci r?) + 
C 1 = 1 


‘2 

k 2 i 


( 20 ) 



13 . 


• * 

where 3. and y. are the time rates of change of the variational 

variables corresponding to 3^ and This type of damping on the 

end bodies could be accomplished by employing three orthogonal mag- 

2 

netometers together with appropriate electronics. It is also 
assumed that the cable will provide viscous damping proportional to 
the cable rate of extension. 

III. B. Examples of Terms Expressed in the Nine Variables 

The angular velocity of a body is expressible in terms of the 
angles used to describe the transformation matrix and the time rate 
of change of these angles. From the transformation matrix, equation 
(la), 

®*B/A " (flC0SB3C0S32 + ^2 s i n33)b-j 

• > • * 

+ (-3]Sin3 3 cosg2 + 3 2 cose 3^2 + (3ySin3 3 + 3 3 )b 3 (21) 

It is noted that, in accordance with the assumed one-two-three 
sequence of rotations, equations (la), that corresponds to a 
rotation about the A-j axis, gg corresponds to a rotation about the 
displaced ^ axis, and £ 3 is associated with a rotation about the 
Bg axis. 

The expression for u", the relative translational velocity vec- 
tor, can be reduced for planar motion to the following: 

A 

u = [i - + 52 ) sin (B 3 - ep -p 2 (y 3 + sin (y 3 - e 1 )]d 1 

+ [(e-j + a)l + ppSj +n)cos(63 ' e i ) + P2^’ y 3 + 52 )cos(y 3 - e-|)]d 2 ^ 22 \ 
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i 

i 

-6 

while Stabekis has previously developed for the same quantity: 
u = [i - p-j (e'j + i -j )si n4> -j - pgCo-j + & + ^^^^l 

a. 

+ [(e-|+ &)*,+ p-j(o-{ + J2 + ii)cos^ + P2^l + ^ + ^2 ^ cos ^2^2 

By comparing the two, we can relate to similarly i> ^ 

to y 3 

3 3 = e-j + <j>f Y3 = e l + ^2 

for planar motion. The angles ^ and <j> 2 describe the orientation 
between the cable line and the principal axis of each end body which 
is aligned with the attachment arm vector (Fig. 3). 


III. C. Lagran ge's General Equation.._a_n d .t he . P r o c edure, for, 
' Developing the Equations of Motion 

With all expressions, kinetic, potential energies, etc., in 

terms of the nine generalized coordinates, Lagrange's equations, 

_d ( 9T \ _ 3(T-V) = _ 3f_ (25) 

dt 3q j 3q j 3q . 

yield the nine equations of motion for this space $tation--counter— 
weight system. 

The expansion of equation (25) can be simplified by using index 
notation. For example, 


a(T-V) _ l i v u. 3u i 

i=r 1 


+ 1 






Bt W B- 


3w B-j 

8q 


+ I r w r 
c i c i 


3w C-j 

"3qT 
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’Sqj 


V* 


S-Yj 

3qj 


)- 


%t 


( 26 ) 


The total expression for T-V need not be expanded before differ- 
entiation. Similar expressions were used in evaluating 3 T/ 3 qj. 

III. D. The Equations of Motion 

The nine equations of motion were derived with the approximation 
that sin(qp) 5 qp and cos(qp) s 1 where q p is any one of &■) , 32, ©2- 
Y*j or Y2 (angles out of the orbit plane). The approximation was made 
after the final differentiation and all terms of degree higher than 
two in q D were considered small. In addition it is assumed that the 

r 

attachment arm vectors r^p, r^/p are direction of the unit 

A ■ A 

vectors, b-j and , respectively; for the more general case where the 
cable attachment point is not located on the B] (C-j) principal axis, 
the equations would have to be appropriately modified. 

One of the nine equations obtained after only these approxima- 
tions appears in the appendix. All of the remaining eight equations 
are more complicated algebraically than the 0^ equation referred to 
here. 


III. E. Linearization of the Equations of Motion 

In a stability analysis concerned with motion about an equilib- 
rium state, variables are used which measure the deviation from the 
equilibrium motion. The following definitions were used in this 
respect. 
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j 0-j = 0 n t + X <$ = 

e 3 = V + “l r 3 = V + a 2 W 

where e n t is the equilibrium value of G], 33, and Y3 for any time, t, 
and is the equilibrium cable length to be determined. x s , <*£, 

and 6 are the new variational coordinates corresponding to 9], 63, Y3 
and £, the original variables. The original out-of-plane angles are 
zero at equilibrium and accordingly serve as variational coordinates. 
The variational form of the 9] equation in which terms of degree 
higher than two in the variational variables have been neglected is 
shown in the appendix. 

After examining the variational form of the equations, the 
linear equations are obtained by making use of trigonometric identi- 
ties and neglecting all second and higher degree terms in the varia- 
tional variables. The approximations 

sin(33 " e l) = sin(a-j- x) ~ orj- x and 

cos(P3 ~ 9-j) « cos(a-j - x) - 1 , etc., also yield 

linear terms. The linear equations upon which the subsequent stabil- 
ity analysis is based are shown as follows: 

X~ X + P 1 « 1 + P 2 « 2 + 2(® n + «)& + (P-J + P 2 )(®n + 

1 * 2 . * p 

-p-j(e n + ft) a-j - p2(®n + ft) 


a 2 - 0 


(28a) 
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6 ~ 


6-0 + p-| + p 2 ) ( ®n + “2(6 n + a)x - 2 p i ( e n + n)“i 


i . 


-2p 2 ( + “ (^n * 6 + M 6 " 6 o) + k 2 5 55 0 ( 2Sb ) 

I rt 1 «i If** 1 ^ ^ » • 

a-j — (p-j + Ig^) a l + P 1 p 2 a 2 + + 3 °! + 2 p l( 6 n * 


-P 1 (e n + fi) 2 x+P^(l + P 2 )(e n + fi) 2 « 1 +c B 3 o‘ 1 -P 1 P 2 { 0 n +«) 2 “ 2 = 0 < 28c ) 

V( p 2 2 + Ic 3 )°2 + p l P 2 a l + P 2 X + k C 3 “2 + 2P 2 ( 0 n + «)* 


-P2(6n+fi) 2 X+P2(l +p i)(Sn + ^) 2 “2 + cj. “ 2 - p -] p 2 (® n +JI ) 2a i = 0 ( 28d ) 


h“ 


h- 


*| j 

1 Hb 2 + ^ h +[{I B 2 + I B 1 )n - I B 3 t ft n + ’ n ^ 6 l + V2 

- 2 I B 3 + V n 2 3*2 = 0 ^ 

' 1 ^B 2 + + + ' + I B 1 ^3 h + k B 1 ®V 


2 D 1 


+[<*, + -I < r B n - 2I B, + Z bX 3 = 0 < Z9b > 


Bi B 3 B 2 


e 2 — e 2 + (1 + p -|+ p 2 )(0 n + a) e 2 - 0 


( 29 c) 
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t i 


Y Z ~ I dc 2 + V " 2 +[(I C 2 + V'° - V 9 " + a)] *i + V 2 


+ C c C + T C 3 V " 2 (I C 1 ' 2I C 3 + l C^ ^2 " ° 


Y r- \ (I C 2 + ' Y 1 + [1 c 3 < S n + - (I C 2 + l 'c^2 + k C! Y l 


■[Cc + Ic 3 - 2I c 3 + h 2 >° 2 3*1 =0 


(29d) 


(29e) 



The linear equations separate into four parts. The in-plane 
equations (28a)-(28d) are completely uncoupled from the out-of-plane 
equations. The out-of-plane equations are made up of three separate 
parts: (1), the 82 , g] equations (29a) and (29b) which are dynami- 

cally coupled, (2), the 82 equation, (29c), and (3), the Y 2 > Y] equa- 
tions, (29d) and (2Se) also dynamically coupled. The $1 equations 
apply to body 1 in the same way that the y 2 , y-j equations apply to 
body 2. The equation, however, separates completely from the other 
equations and indicates simple harmonic motion. Although asymptotic 
stability of the system clearly will not occur for the case in which 
the cable has an out-of-plane perturbation, mission requirements could 
perhaps still be accomplished for small disturbances. From the defi- 
nition of the angle 02 , the system would achieve an equilibrium 
•motion in a plane slightly inclined with respect to the original 
plane of motion, but with the same equilibrium length and nominal 
spin-rate. For such a situation, all stability criteria previously 

developed by considering the original equilibrium motion would apply 

9 

to the new equilibrium motion. 

At equilibrium, the values of the variational coordinates are 
zero. An examination of the l equation at equilibrium yields the 
equilibrium condition 
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With 6 0 = ■ and p- = p^/Z & for i =1,2, equation (30) becomes 

k^,, - l 0 ) = u(£ e + pi + P2)(«n + n > 2 (31) 

which states that the centrifugal force due to the inertial rotation 

of the system is balanced by the tension in the cable. Solving 

equation (31) for £ e reveals that 

p _ hA) + v(p 7 + P2^®n + (321 

e h - 

The condition k-j ^ y(6 n + ft) 2 implies that l c is either negative or 
infinite. For realistic values of it is necessary that k-j be 
greater than y(e n + ft) 2 since if this condition is not satisfied the 
cable will not provide sufficient tension to balance the centrifugal 
force of the spin. 

The fact that the in-plane equations all separate from the out- 

of-plane equation allows a means of comparing the in-plane equations 

6 

with those in Ref. 6. The analysis of Stabekis was confined to the 
orbit plane. From equations (24) and the definition of the varia- 
tional variables, equations (27), one can relate the variational 

coordinates of Ref. 6 with those used in the present analysis: 

a-j = X + <j>^ = X + <f>2 (33) 

The in-plane linear equations can thus be written in terms of the 
variables defined by Stabekis. 

After using equation (33) and manipulating the in-plane equa- 
tions, equations (28a) - (28d), it was found that equation (16) of 
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I 

! 

Ref. 6 did not include the linear terms yp^-j f^+p^+p^) > and 

yp 2^2^ +p l +p 2^ while in equation (17), the terms -2p-j (e n +n)^-j , and 
**2p2(B n +ft)<f>2 v;ere missing. The equations of motion for the system 
used by Ref. 6 were completely rederived using the variables defined 
therein and confirmed the indication that these terms were missing. 
When converted to nondimensional form, these four terms all contain 
the coefficient p -[/£ c which is much less than one for the examples 
considered (attachment arm lengths much shorter than cable equili- 
brium length). Thus the effect of neglecting these terms on the 

6 

numerical results previously reported would be expected to be small. 
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IV. STABILITY ANALYSIS 

IV. A. The General Stability Criteria 

■v The procedure used in analyzing the stability of the general 
system was to first obtain the characteristic polynomial. The coeffi- 
cients of this polynomial can be considered in conjunction with the 
Routh-Hurwi tz necessary sufficient conditions for the roots to have 
negative real parts, in order to provide conditions for stability. 

The in-plane and out-of-plane criteria can be obtained independently 
of each other because of the decoupling of the equations. 

The first step in obtaining the characteristic polynomial is to 
nondimensional i ze time as follows: 

t = (e n + a)t (34) 

This simplies the coefficients of the variables in the linear equa- 
tions. Differentiation with respect to time is now replaced by differ 
entiation with respect to t and primes on the variables replace dots. 
(Note that primes above the various coefficients represent nondimen- 
sional ization of physical parameters and should not be confused with 
differentiation. ) The nondimensional linear equations are shown 
below: 

H I II I II Il« I I 

X— x + P-jCt] + P2 a 2 + + (p-| + P2)x-p-} a ]-P2 a 2 = 0 (35a } 

ti i ii it. H it i 

<$-- <$ - 2x - 2 Pf ct] - 2 p 2 a 2 + + ^2° = 0 

1 2 i - it i i ti t ii u ii 

“l" + + p l p 2 a 2 + p l x + * < B 3 C1 1 + 2p 1 5 


(35b) 
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-p'x + p’(l+ P 2 ) a 1 +-c B V P l P 2 a 2 = 0 


a 2 - 


32- 


(p‘ 2 + l c y 2 + p ] p 2 a l + P 2 X ” + k C 3 a 2 + 2 p 2 6 

, i i " * ’ 

-P 2 X + P 2 (l + p l) a 2 + C c 3 a 2 ' p l p 2 a l - 0 

\ (Ib + i b ’)B2 + [( j b 2 + - i b 3 3bi + k « z 6 2 


+[CB 2 + I B 3 r - 1 Pb! + i B 2 ) 1 ’ 2 ] 62 - o 
B r - £ (Ig 2 + + [Ib 3 - ( j b 2 + I Bi) r ^ s 2 + k B-j 


+ [c B 1 + I B 3 r - \ (I B 1 + I B 2 )r<: ] 6 1 = 0 

It * 1 

e 2 — 02 + (1 + p l + 9 V e 2 “ 0 


T 2 - \ dc 2 + hy)y Z + Cdc 2 + ] C! )r - ^3^1 + k C 2 Y 2 
+[Cc 2 + Ic 3 r " \ (I Ct + ]y 2 = 0 


It it I 


Y 2 + Vi 


+Ecc-j + ic 3 r - \ (ic^ lc 2 )r 2 hi = 0 


(35c) 

(35d) 

(36a) 

(36b) 

(36c) 

(36d) 

(36e) 



24 . 


Considering first the in-plane equations, equations (35a) - 
(35d) , the substitutions q p = A r e xt , where q r is one of the in-plane 
variables (r = 1,...,4), yields the following equation necessary for 

a nontrivial solution for the A r ‘s: 


X 2 + R 

2 A 

P* (x 2 -l ) 

p'(X 2 -l) 

11 


1 

2 

-2X 

0 11 ,11 

x^+k x+k 

i 

-2p'x 

2 1 

“2p^X 

2 

(x z -l) 

2p^ X 

P 33 x2+k B3 X+R 33 

» « . ? 
V 2 (x 


2p’x 

p 1 p 2 (x2_1 ) 

P 44 X +k c x ' 

I (37) is 

an eighth order polynomial in 

X where 


R n ~ p i + p 2 


R 33 = C B 3 + p l° + p 2 ) 

It I « 

R 44 = c C 3 + p 20 + p l) 
'2 1 

P 33 =P 1 +I B 3 


= 0 
(37) 


I O » 


P 44 = p 2 + X C 3 

Expanding equation (37) and collecting terms of like powers in X, the 
characteristic equation for the in-plane motion results, 
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a Q X 8 + a l X? + a 2 x6 + a 3 x5 + a 4 x4 + a 5 x3 + a 6 X? + a 7 X+ a 8 = 0 

The coefficients appearing in equation (38) are related to the system 
parameters as follows: 



a l = Ib 3 l c/z + I C 3 k B 3 + l %/c 3 

a 2 = (1 + p ! + P2KP33P44 “ P 1 2p 2 2 ) + ( R 33 + P 1 2 ^C 3 + k B 3 k C 3 

+ ( R 44 +p 2 ^B 3 + k 2^ k B 3 I C 3 + k C 3 ! B 3 ^ 

+ (ki - 1) i b 3 I C 3 + 4I B 3 ! C 3 

a 3 = kg 3 [(P] + p 2 ) p 44 + ( k l + ^) r 44 + R 44 - 2P 2 2 ~ p 2 2k l3 


+ k C 3 « P l + P 2^ P 33 + ^ k l + 4 ^ P 33 + R 33 


2P^ 2 - P' 2 kp 


+ ^2tlB 3 ( R 44 + 2 p 2 2 ) + i C 3 C r 33 + 2p l 2 ) 


+ (p] + P2 )( r 33 p 44 " P V P 2 2 )] + k 2 k B 3 k C 3 
a 4 = O + p; + P 2>tR 33 P 44 + P 33 R 44 + kB 3 k C 3 + 


+ M k B3 P 44 + k C3 P 33) + k l( p 33 p 44 ’ P 1 2 p 2 2 ) ] + ( R 33+P] *) ( 844 +^) 
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+ k 2[ k Bo( R 44 +p 2^ +k Co( R 33 +p 1^ + ( k l~'* ) [C R 44 +P 2^^B' 


+ kB 3 kC 3 + (^33 + ^1 ^) IC 3 D-^2( kB 3 IC 3 + k C 3 lB 3 ) 


In I II 11 I II I 


- + 4 C( R 33 +p l 2 ( 2 +p i +p 2)) I C 3 + k B 3 k C 3 


+ (R44 + p^ 2 ) 1 ^ + 4p 2 Z ^ + P 1 + p 2^B 3 


= ^b 3 ^ P- | + p 2^44 + ( k l + 4 ) R 44 + ^ p 2 2 + 2p 2 2k l 


+ (p^ + ^ 2 ) ^-1 ^*44 + 4 (p*| + P2^ p 2^3 


+ fcc 3 [ (p l + p 2 )R 33 + (k l + 4)R 33 + 7p i 2 + 2p l 2k l 

I I II I lip 

+ (p-| + P 2 ) k-j P 33 + 4(p] + P 2 )P] ] 


+ k 2 [{P] + p 2) r 33 p 44 + C p ] + p 2) r 44 p 33 + 2 p 2 2r 33 


ip 1 p 1 i p r 1 o t 9 1 ,l 9'9 

+ 2P^ R44 - p 2 iBo - P 1 ^Co + 3P 1 p 2 + 2 ( p l + p 2) p 1 p 2 


+ R 33 R 44.] + C p ] + p 2^ k 2 k B 3 k C 3 


1 o 1 O 11 II 


= (1 + p 1 + p 2 )[R 33 R 44 - P T 2p 2 2 + k 2^ k B 3 R 44 + k C 3 R 33^ 


11 11 11 1 p 1 p 

+ k] ( R 33 p 44 + P 33 R 44 + k B 3 k C 3 + 2p 1 p 2 
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+ (k^ - 1 KR33 + P] 2 )(R44 + P2 ? ) " k 2 £ k B -( R 44 + p 2 ) 


+ k C 3 < R 33 + p l 2 ^ - k l E< R 44 + p 2" }I B 3 + k B 3 k C 3 


101 hm 


+ (R 33 + P ! 2)I C 3 ] +4 (R 33 + P i 2)(R 44 + P 2 2) 


1 2 1 1 1 2 1 2 * ‘ * 2 

+4p] (] + P-j + P2)( r 44 + p 2 ) + 4 p 2 (1 + p ] + P 2^ R 33 + P 1 ' 


I . H 1*3 II, 


a 7 = k B 3 C(P l + P 2 )k l R 44 ' P 2 k l 3 


11 1 • it 1 2 " 

+ k CoCf p l + p 2) k 1 R 33 “ P 1 M 


n 1 1 


+ 4[( P ; + P 2 )(R33 R 44 - p i 2p 2 2 ) - P i 2R 44 " p 2 2r 33 ' 2p l 2p 2^ 


H _ . I > . M " « 1 . _ ‘ 1 \ / “ 1 “ M 

a 8 = k l C( p l + P 2 )c B 3 c C 3 + p l p 2 (1 + P 1 + P 2 )(C B 3 C C 3 )] 
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Because of the complexity of the coefficients in equation (38), the 

necessary and sufficient Routh-Hurwi tz criteria were not developed 

for the general case; nevertheless, equation (38) was used to obtain 

the roots numerically and to obtain graphs corresponding to those 

6 

acquired by Stabekis in his optimization procedure. In addition, 
analytic stability criteria were developed for special cases, and 
for the genera] case, by considering the signs of selected coeffi- 
cients in the characteristic equation. 

The necessary condition for stability is that all of the coeffi- 
cients in the characteristic equation have the same sign and be non- 
zero. By inspection of the coefficients in equation (38), it is 
seen that in-plane stability is insured if at least one of the fol- 
lowing forms of damping is present: cable damping (kp, 0 r rotary 

damping in at least one end body about an axis nominally perpendicu- 
lar to the. plane of rotation (kg^ or k^). Under this condition all 
of the odd coefficients will be positive non-zero. From considera- 


tion of a^, a restoring torque capability must be present on at least 
one of the end bodies about an axis perpendicular to the nominal 

it 

plane of rotation— ei ther c B or c c ^ > 0. Also for a g > 0, ^ must be 


greater than zero. Thus it must be true that k-j be greater than 
u(G n +Q)^ j n order to allow the possibility of stabil i ty--in agreement 


with the results obtained from the equilibrium condition, (viz. dis- 


cussion in connection with equations (30) - (32)). 
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/ 

/ 

The out-of-plane equations (only the equations for one end body 
need be considered, e.g. equations (36a) and (36b), since the equa- 
tions for the second end body are analogous to the equations for the 
first end body) yield the following characteristic determinant: 

=0 (39) 

where 

n 11 = <n 22 = \ (Ig 2 + 1^) 

C 12 = ~ C 21 = I B 3 ' Ob 2 + 

K 11 = C B-| + ! B 3 r " ?Ob 2 + T B ^ r2 

K 22 = c B 2 + T B 3 r ‘ ¥ l &2 + 

The characteristic equation for this fourth order system is: 
m ll x4 + + k B 2 ^ 3 + Cm l 1 ( K n + K 22^ + + c l2^ x2 

+ ( k B/22 + k B 2 K n^ x + K 11 K 22 = 0 <«» 

which can be written, 

b 0 A^ + b-j + ^> 2 ^ + ^A + b^ = 0 



(40a) 
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The Routh-Hurwi tz necessary and sufficient criteria are that for 
b 0 > 0, it is true that 

A. bj > 0 

B. b-j b 0 

= b*ib 2 - b 3 b 0 > 0 

b 3 b 2 

C. b*j b 0 0 

2 

b 3 b 2 b] = b-| ( b 3 ** b*jb^) - b^b^ > 0 

0 b 4 b 3 

b-j b 0 0 0 

b 3 b 2 b-, b 0 

>0 

o b 4 b 3 b 2 

0 0 0 b 4 

The fourth condition is satisfied if the third condition holds and, 
in addition, b 4 > 0. when the b * s are written in terms of the para- 
meters, the four criteria for roots with negative real parts are 
obtained as shown below: 

A. ^B] + k B 2 > 0 

B - m n k B 1 c B 1 + m n k B 2 c B 2 

+ ( k By + k B 2 )( k B 1 k B 2 + Ob 3 - j (Ib 2 + + m fl r2 ) > 
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m llO<B 


ke 2^ lB 3 


Ob 2 + 


B,) r )r 


II ii t ii ii . o r , » O II , Ho h 

C. m ll k B 1 k B 2 ( C B 1 ' c B 2 ^ + C k B 2 C B-, + k B^ C B 2 

, . II , Ii .0,1 1 I l v x , H , H . It II i-tr. Ii , H 

+ (VV ( V? ( VV r) b iVV Cb 2 )][ V b 2 


+ 0b 3 - (Ib 2 + ^ ] > 0 


D. [c" i+ (I' - 1 (I‘ 2+ I^Dr ][cb 2+ (I' 3 4 (le/I^Dr] > 0 


From condition A, there must be rotational damping on each of the end 
bodies about at least one of the principal axes which, at equilibrium, 
will lie in the plane of rotation. {Note that rotational damping must 
be present on both end bodies since similar criteria may be developed 
for the second end body). Furthermore from condition B if 

E. eg. i Ig^r - \ (Ig 2 + I B ’)r 2 , for both i =1,2 

asymptotic stability of the out of plane motion is assured if condi- 
tion A is also satisfied. It is also seen that the satisfaction of 
E guarantees condition D. 


IV. B. Special Cases of the Linear In-Plane Equations 

Assumptions on the physical model can reduce the complexity of 
the in-plane stability analysis. The cases of- identical end bodies; 
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where the cable is attached to the c.m. of both end bodies; point 
mass end masses; and the case in which the system moves in orbit 
with a constant inertial orientation are treated. 

The Case of Completely Identical End Masses 
For this case, it was assumed that the space station-counter- 
weight system was completely identical, that is, that p-j = P 2 > 
m-j = m 2 > snd for i — 1,2,3, Ig, — ^ q . » ^B* ~ ^C* ~ ^C • * 

The characteristic equation separates into two factors for this case 
as shown below: 

f ( a) • g(x) =0 

where 

i 1 nil 

f(x) = Ib 3 X + kg 3 X + R + P-| 2 ; R = Cg 3 + p-,(l + Pj) (41) 
and, 

g(x) = IB3A 6 + (^83 + ^2^3)^ + f c B3 + p l(" 1 + 2P ] 2 ) 2 

+ + (k-j + 2 p-j+ 4) 

+ ^(cg^ + P|(l + 2 p-| 2 ) )+kg 3 (k-j+2p-{+4)+2P-jk2lg 3 ]A 3 

+ [2p-j (1 + 2p^ ) + (k-j + 2 p-j + 4)(c B ^ + p l^ + 2 p-j ) ) 

1 11 11 • 11 1 2 

+ 2P^k 2 kg 3 + ^P-jk-jIg^lx 
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+ [2p^ 2 (1+2p 1 2 ) k^P^ ( c b 3 +p i 0 +2p ] 2 ) ) +2p i k-jkg^A 


to I O 9 tl I It !l I • ? 2 

+ 2 p ] z (1 + 2 p ^) k 1 + + <^0 + zp 3 J ) 


(42) 


The quadratic factor, f(x), indicates roots given by 

X = ± _L [k^ - 4 Ib 3 (R + P1 2 )] V2 

2I B 3 3 


(43) 


where R = Cg^ + p-j ( I+p-j ) - In-plane instability is associated with 
this mode for kg 3 < 0 or s - p \( 1 + 2 p \) . However, it should be 


recalled from consideration of equation (38) that either c B ^ or 

> 0 is required for stability of the general system--i .e. 5 a stronger 
criterion on the restoring constant than is apparent from equation 
(43). The .fact that the missing terms of Ref. 6 are not involved in 
this factor makes the results obtained above analogous to those 
obtained therein. 


The Case of Zero Attachment Arms 

i « 

When p-| = P 2 = P-j = P 2 = ^ the ca ^ e 1S attached at the c.m.’s 
of both end bodies and the in-plane characteristic equation separates 
into the factors:' 

l n n it 

i B 3 x + k B 3 x + c B 3 


( 44 ) 
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‘c/ + k C 3 X + C c 3 


and x 2 (x 2 + + k-j + 4) (46) 

The repeated zero root resulting from expression (46) is indicative 
that in-plane asymptotic stability never occurs for this case (see 
the appendix). But aside from this, stability is indicated in the 
other modes for c^ > 0 and > 0; c^ > 0 and > 0. These 
results indicate that rotational damping and restoring effects must 
be present on both end bodies, as well as cable damping and restoring 
forces. This is a stronger criterion for in-plane stability than 
that deduced earlier in connection with the sign of the coefficients 
in equation (38). 

For an actual situation where the attachment arms are much 
shorter than the cable length p. = p./£ £ * 0, the results obtained 

here would have important implications on stability. 


The Case of Point End Masses 

i i 

For this case, pj = p£ = 0 and there is no rotational end body 
motion so that only the l and 8-j equations remain. The in-plane 

characteristic equation will contain only those terms appearing in 
expression (46) and the repeated zero roots still will occur indicat- 
ing the presence of a first integral of the motion. 
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The Case in which o n - -n 

Here, the system spin rate is equal and opposite to the orbital 
spin rate so that the system moves around the orbit with a constant 
inertial orientation. The 0 ] equation involves only the acceleration 
x , and a 2 . When the 0-| equation is used to eliminate X from the 
a and a 2 equations, the ^ and <* 2 equations resulting are exactly 
those obtained for p’ = P 2 = 0. The roots indicating the cable 
motion are obtained from 

X 2 + k^x + k.j = 0 (47) 

which indicates asymptotic stability in this mode for k 2 > 0 and 

> o. A more conclusive discussion of this limited case would have 
to include the effects of gravity-gradient torques, which become a 
more important perturbation on the system motion for lower spin rates 


IV. C. Special Cases of the Out-of-Plane Linear Equations 

The out-of-plane equations, yielding a fourth degree characteris- 
tic equation, allowed the determination of the Routh-Hurwi tz conditions 
for the general system. The next most simple case is that for which 
the 32 anc5 ^1 equations uncouple (C-j 2 = Cg] - 0 in equation (39) ). 

In this instance, 


^83 

T B 2 + r Bl 


r 


(48) 


The 3 2 anc ^ equations become two second order equations from which 
it is clear that asyptotic stability occurs in these modes for 
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and c 


B. 


i - 


_ Jl f t * 
2 ■ B. 


+ Ii )r‘ 


(49) 


with kg^ and kg^ both positive. 

For a realistic system where the spin rate is much greater than 
the orbital rate, r « 1. Therefore equation (48) would be satisfied 
only by bodies having Ig^ « Ig^ + ^B] v/ ben the enc ^ bodies have their 
"long axis" perpendicular to the nominal plane of rotation, (i.e. for 
long, slender pencil -shaped bodies). 


V. FIRST ORDER GRAVITY-GRADIENT EFFECTS 


The following symbols will be used in this section: 

^ unit vector in the direction of the system local vertical vector 

G the gravitational constant 

1^ the inertia dyad of body 1 with respect to the B system 

1^ the inertia dyad of body 2 with respect to the C system 

M e the mass of the Earth 



unit vector corresponding to the direction of increasing qj 
generalized force associated with the qj coordinate 


R 

\ 


the constant orbital radius 

the torque due to gravity-gradient effects on the first end body 
(the body in which the B system is fixed) about its center of mass 


T(; g the torque due to gravity-gradient effects on the second end body 

(the body in which the C system is fixed) about its center of mass 
7 

Robe used the following expression for the gravity-gradient torque about 
the mass center of the first end body of a tether connected two body gravita- 
tionally stabilized system: 


Tb 


9 


3GM e 

2R 3 


(§1 x i B • a-|) 


Similarly, for the second end body, 


( 50 ) 
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It should be noted that these expressions are first order approxima- 
7 

tions and would not be valid for the case of very long separation distances 
between the two end bodies. 

Under the same assumptions used to obtain the equations of motion 
(equations (28) and (29)), TBg and Tc g can be written in terms of the nine 
independent coordinates of this system. If N is the total gravitational 
torque on the system, then it can be shown that 

N = TB g + Tc g = T 2 a 2 + t 3^3‘ ( 52 ) 

Assuming all second degree and higher terms in the out-of-plane angular 
coordinates are small, the expressions for T 2 and T 3 can be developed as 
follows: 

T 2 « ^§[-(I Bl - Ib 2 )M si " 2(e„t + ai) 

+ OB] " 2I B 3 + i B 2 ^ 3 2 + ( I B’( " i B 2 ^2 cos 2 (®n t + a l) 

" ^ I C 1 “ V Y 1 Sln 2 ^ nt + + ^ C 1 " 2I °3 + Ic 2^ 2 

+ “ V T 2 cos 2 ( ®n t + a 2^ (53) 

and 

3GMe r » 

t 3 = _ 2p^ I B 1 ' I B 2 ^ sin 2 ^ e n t + “P 


+ Oc-| - 2 ( 6 n t + a 2)] 


(54) 
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The following transformation is applied, according to the principle 
of virtual work, to convert the gravitational torque into generalized 
forces: Q q = n q ■ N". Lagrange's equations then become: 

5 j 


A (Ah a (T-V) = n _ aF 
dt V 3 <io " qj 3 «j 


(55) 


In equation (55), V is the potential energy without including gravity- 
gradient effects. 

To a first order approximation, the gravitational force on body 1, 
Fgg» is given by Ref. 7 as: 
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Similarly, the gravitational force on body 2 can be written: 



GM £ m2 

T 


CO - 


r c* h 

R 




(56b) 


where Fg is the vector from the system center of mass to the center of 

mass of body 1 and Fq is the vector from the center of mass of body 2 to 

the system center of mass. Equations (56a) and (56b) can be used to show 

that the total gravitational force on the space station-counterweight 

system is, to first order, not a function of the relative position vectors, 

7 

rg and r^ . Therefore there will be no first order effect of these forces 
on the equations of motion (equations (28) and (29)). 

The first order generalized forces due to gravity-gradient effects will 
now be evaluated. The n q ^ for the angular coordinates are developed as 
(viz. equations (1) and (2)): 
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n S 1 = a i 


"Sz = 

cosB-j^ + sinB-j^ 


"e 3 = 

s i n02 a l “ cos$2sinS'|a2 + cosB2CO$B]a3 

' (57a) 

A 

% = 

h 


‘ A 

n v = 
2 

cosV^a 2 + sinY^a^ 


A 

n y = 

’3 

sinY2a-j - cosY2sinY-ja2 + cosY2CosY-ja3 

(57b) 
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Thus the generalized forces, Q q . = fi q • N, under the assumption of small 

amplitude displacements in the out-of-plane coordinates, can be expressed 
as follows: 

Qg, = a i • n = o 

% 2 ‘ T 2 + ®1 t 3 

Q 63 ' t 3 (58a) 

Qy, = 0 . . 

Qy 2 ' T 2 + Y lT3 

Q Y 3 ' T 3 (58b) 
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-cose-|T 2 


(58c) 


It should be noted that = Q y ^ =0 since according to equation (52), 

IT has no a-j component. Similarly, Q 0 ^ and Qq^ are also exact expres- 
sions within the assumptions previously stated. 

Equations (28) and (29) are in the dimensions: rad/sec 2 because they 

2 

were obtained from Lagrange's general equations after division by \il z . 

(an exception is the t equation (equation (28b)) which was obtained from 

the general form of Lagrange's equation by division by y). For dimensional 

2 

consistancy, equations (53) and (54) must also be divided by y£ e . The 
resulting expressions are: 


» 3GMp « » * 

t 2 = HIb] - I B 2 ) & l s1n 2 < 9 n t + a l) 

+ Ub*, - 2Ib 3 + h z )$Z + (Ib-j “ i B 2 ^2 C0S 2( e n t + a l) 

- 0(/| “ I C 2 ^l sin 2 ( 9 n t + a 2) + (I Cl " 2I C 3 + ! C 2 ^2 

+ (I c - I(O y 2 c0S 2(0 n t + a 2 )] (59) 

I 

T 3 " COb, - I B 2 )sin 2 ( 8 n t + °1> 

+ ^ C 1 " V s1n 2 ^n t + a 2^ (60) 



42 . 


which are in the units of rad/sec^. 

Thus to a first order approximation, the linear equations of motion 
with gravity- gradient effects are: 


X— 


x + P-J ot-j + P 2 c*2 + 2(0 n + 


i i* 2 i* 2 i * 2 

+ (p-j + P 2 )(0 n + x “ p l ( 0 n + a) <*1 - p 2( 0 n + «) a 2 


= t 3 


(61a) 
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(61b) 


i / i .. 


“i 


( p i + + p i p 2<%2 + PiX + kg 2 a l 


i i 


■ t 2 , 2 

+2p*j ( e n + £})<$ - Pi (e n + ft) x + p ^ ( 1 + p 2 ) ( ® n + ^ a l 


Cg ot-j - p i ( e n a 2 1*3 

3 


(61c) 


1 2 


do — 


{p 2 p l p 2 a 1 P 2 X+ k c oi 2 


1 • . 1 . 2 * 1 * 2 

+ 2P2(Q n + °) 5 " p 2( 0 n + «) ^ + p 2 (1 + P])(e n + n) c <2 


+ 00^2 - Pl p 2( 0 n + ^r a l “ t 3 


( 61 d ) 



43 . 


$2- 


\ (Ib 2 + + t(h 2 + %)* 


i i * 


-I B3 ( 0 n + n)]Bi + k Bz $ 2 + Ccb 2 + I B 3 e n fi 


B l" 


2 Ubi ' 2I B 3 + I B 2 )fi ]e 2 = T 2 + 6 1 T 3 

\ Ub 2 + + [I B 3 (9 n + °) 

■(Ig 2 + ^-j) ^2 + k Bi e l + ^ C B-) + 1 B3 9 n Q 
' ^B] " 2I B 3 + r B 2 ^ ^ e l = 0 


(62a) 


(62b) 


62 + (1 + p i + p 2) ( ® n + n) e 2 - - T2 cosQi ( 62 c) 
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It can be seen by examining equations (59) and (60) that the first 
order linear equations with gravity-gradient effects (equations (61) and 
(62)) now involve periodic coefficients with frequency at twice the spin 
rate. Also, each of the in-plane equations (equations (61)) except the 
t ( 6 ) equation, now contain forcing terms of constant amplitude on the 
right side with frequency 2 e n . 

Certain conclusions can be drawn from the expressions for T 2 and 
T 3 (equations (59) and (60)). First, as noted previously, the in-plane 
motion is forced by terms of frequency equal to 2e n . Second, it is 
apparent that gravity-gradient effects become more pronounced for 
small R. Also these effects are increased as either 
Ubt - !b 2 I or l 1 ^ - ! C 2 I 
or both are increased. 

A rigorous stability analysis of a linear system with periodic coeffi- 

. 10 

cients could be made using the Floquet theory. For a complex system, the 
application of the Floquet analysis would necessitate extensive computer 
simulation to examine the moduli of the eigenvalues associated with an 
augmented matrix and evaluated over a wide range of system parameters. 
Although the Floquet theory was not applied to this study, the effect of 
gravity-gradient torques on the system was considered numerically for 
selected steady state responses as well as transient responses. Resonance 
due to gravity-gradient effects is shown to exist for certain special 
cases, easily identified, and these results are presented in Section VI. C. 
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VI. NUMERICAL ANALYSIS 

When the least damped mode of a system is examined, the worst 
possible response of the system is examined. All of the other modes 
decay at a faster rate than the least damped one. If it can be 
determined for which set of system parameters the least damped mode 
decays fastest, then the system damping can be optimized. The roots 
of the system characteristic equation can be calculated numerically 
for a specific set of system parameters. By incrementing the system 
parameters one at a time, a complete range of physical constants can 
be considered. Then by finding the least damped mode and plotting the 
the time constant of this mode as a function of each parameter, the 
optimum set of system constants can be determined. 

In this section, the numerical results obtained by the procedure 
outlined above are compared with those results given in Ref. 6 for 
the in-plane characteristic equation. Thus only those parameters 
treated in Ref. 6 are considered here. Furthermore, the out-of- 
plane modes are examined for the same range of parameters. Transient 
responses of the linear equations for two different sets of param- 
eters are also presented. 

All computer results were obtained by means of an IBM 1130 
computer system. The roots of the characteristic equations were 
calculated using a Newton-Raphson iteration technique. The transient 
responses were obtained by integrating the linear equations (equations 
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(28) and (29)) employing a Runge-Kutta fourth order method. The in- 
plane optimization program required about fifty minutes for each 
curve of T vs. system parameter which contained 250 points. This can 
be compared with the out-of-plane optimization program which only 
used ten minutes of computer time for the same number of points. For 
the computational time step chosen. (At = 0.25 simulated problem sec- 
onds), the transient responses required thrity seconds for each simu- 

lated problem second. 

In all computations, for both examples considered it was assumed 
^ " c B 2 = cg 3 = cgi = c £ 2 = c £ 3 

and 

tt »• » « 11 “ 

\ = k B 2 = k B 3 = k C 1 = k C 2 = k C 3 . 

VI. A. Identical EndJSodje s. 

Now, the time constant of the least damped mode, T, was calcu- 
lated as a function of c$. and kg. in this first example. The fol- 
lowing parameters remained constant at the values given below: 
fi = 0.055 deg/s k-j s 1000 p/f 

‘ 6 n = 32.0 deg/s = 56 ' 7 P* s/f 

m 1 = m 2 = 600 slugs I Bl = I Cl = 81,000 sl-f 2 

p, = P 2 = 12 f Ig 2 = I C2 = 80,000 sl-f 2 

t 0 = 230 f lB 3 = = 86,400 sl-f 2 



47 . 


In the first study, the rotational damping constant, kg , had the 

constant value 15,500 f*.P * ~L as the. spring constant, Cn , was incre- 

rad 3 

mented. Fig. 4 shows that T as a function of c^ increases slightly 
as cl 3 increases. Also shown is the plot given in Ref. 6 in which a 
rotational restoring constant of 5000 f-p/rad corresponded to the 
minimum value of t for this same set of system constants. The cor- 

ii 

responding out-of-plane graph appears in Fig. 5. Here, as cg^ 

increase, T decreases in a manner which seems asymptotic. It can be 

6 

seen that the minimum time constant achieved by Stabekis for the 
in-plane modes is about two orders of magnitude less than that 
obtained here for the same range of parameters. (It should be noted 
that in Fig. 4 and all subsequent figures, the parameter shown on the 
abscissa is dimensionless. For the nominal system parameters con- 
sidered the conversion to the corresponding dimensional quantity is 
given below the abscissa on each figure.) 

Throughout this section the optimization curves presented in Ref. 
6 were two and sometimes three orders of magnitude smaller than the 
curve obtained in the present work. This difference in results was 
due to the fact that rotational damping and restoring system constants 

of Ref. 6 were defined with respect to angles and angular rates mea- 

• : 1 

sured from the cable line ($s and $s--see Fig. 3) where, on the other 
hand, in this analysis, system constants were defined with respect to 
variational angles and angular rates which include the effect of x and 



48 . 


X according, to equations (33). 

It should be recalled from the stability analysis of the general 
in-plane equation (38), that (or c^) > 0 is necessary for in- 
plane stability. For this reason, time constants associated with zero 
values of rotational restoring constants are not indicated for this 

unstable situation in Figs. 4 (and 14). 

Fig. 6 is a graph of T versus kg 3 in which Cg 3 was held constant 

« 

at 5000 f *p/rad. Near k B3 = 0, the value of T is very large and as 

kg 3 increases T decreases in asymptotic fashion. The out-of-plane 

graph of Fig. 7 has the same characteristics as the in-plane graph but 

with a minimum time constant of about one order of magnitude smaller. 

6 

Included in Fig. 6 is the graph obtained by Stabekis which shows his 
T minimum is reached for a rotational damping constant of approximately 
12,000 f*p‘ s/rad. 

The transient responses of the identical system studied in Figs. 

4 through 7 with c B- . = 5000 f-p/rad and kg,. = 15,500 f-p'S/rad are. 
given in Figs. 8 to 13. The initial conditions used were, zero velo- 
cities in all variables and, 


X = 0 

ag = -0.1 

rad 

If 0.48 f. 

ft] - 0.1 

rad 

ct-|= 0.1 rad 

$2 = 0 



Initially only a response of 100 seconds simulated time was considered. 
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From Fig. 9, the high frequency motion of the cable is seen to 
be greatly damped for the value of cable damping chosen (k 2 « 56.7 

p-s/f). p or parameters chosen, p = = l2f and = 256f, it 

should be noted that p-j and p 2 have the value of about 0.0468 which 
is much less than unity. The results of the special case of zero 
attachment arms indicate that the in-plane equations are weakly 
coupled for very small p-j and p-j . Thus the graph of x shows predom- 
inately a lower frequency motion, and the response of t - shows 
mainly the motion associated with the 'cable, that is, until this 
motion decays and the effects of coupling become of the same order of 
magnitude. This same type of motion is apparent in the graphs of a-j 
and a 2 * The 9 ra Phs of 6] and s 2 indicate damping, but because the out 
of-plane motion is independent of attachment arm lengths, the effects 
seen in the in-plane graphs do not take place. Since the motion of 
e 2 is simple harmonic of frequency (1 + + p 2 ) (0 n + o) rad/s 

(see equation (29c)), the response of the variable e 2 is not shown. 

All of the transient response considered in Figs. 8 - 13 were 
then considered for an extended response time up to 600 seconds in 
order to reveal the damping of the lowest frequency motion. These 
extended responses required about five hours of computer time. The 
results are shown in Figs. 8a, 9a, 10a, 11a, 12a, and 13a. Assuming 
that the lower frequency motion of the in-plane responses is more 
representative of the least damped mode of motion, the time constant 
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of this motion can be calculated by means of the equation 

T = *2 ~ ^ 

ln(y^/y' 2 T 

where y-j is the positive amplitude occurring at t * t-j , for i = 1 »2 
and is measured directly from the transient response of a variational 
coordinate. Using Fig. 8 a, it can be measured that t 2 - ti ~ 367 sec. 
and yj/y 2 - 1.18*, with these values, T “ 2186.37 seconds which is 
approximately 194.68 revolutions for the given inertial spin rate of 
this system (32.055 deg. /sec.). From Fig. 4, at = 5000 f-p/rad 
( c B 3 ~ 0-0008), T, the time constant of the least damped mode, is 
about 197 revolutions. The value of T obtained from the transient 
response in Fig. 8 a is therefore about one per cent different from 
the value predicted for the least damped mode. However, the accuracy 
in measuring T from the transient motion largely depends on the error 
obtained in measuring the ratio y *i /y 2 * ^ rom consideration of the 
other in-plane responses. Figs. 9a, 10a, 11a, it can be seen that the 
measured time constants are within ±4 revolutions of that determined 
from Fig. 8 a. 

The same procedure applied to the 3 -j and 63 motions reveals that 
the time constant of the response (Fig. 12a) is about 4.35 revolu- 
tions and the time constant of the s 2 response (Fig. 13a) is about 
3.79 revolutions. Since these two motions are more highly coupled 
than the in-plane motions, the transient response do not indicate 
the least damped mode to the same degree as the in-plane transient 
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responses. Nevertheless, the time constants calculated are less 
than the time constant of the least damped mode (t = 7 rev.) obtained 
from Fig. 5 at c^ * 0.0008. 

Thus from the measured time constants based on the transient 
responses of both in-plane and out-of-plane coordinates, it can b_ 
observed that these results are consistent with those predicted by an 
examination of the roots of the system characteristic equation. 

VI. B. Unidentical End Bodies 

For this second example, the nonvarying system constants were 
selected as: 

n = 0.055 deg/s Ib^ = 130,000 sl'f 2 

e n = 32.0 deg/s Igg " 100,000 sl-f 2 

2 

m-j » 770 slugs I B 3 “ 173,250 sl'f 

^2 = 430 slugs IC] ~ 30,000 sl*f 

lo = 230 f IC2 = 20,000 sl'f 2 

P] = 15 f IC3 = 34,830 sl-f 2 

P2 = 9 f 

For this case in which the end bodies are dissimilar, the varying 

” » M f “ j ■ 

parameters were c^., kg^, k-j , and 

In' the first study presented in this section the values 
kg.. = 10,000 f.p. s/rad 
k-j = 1 ,000 p/f 

k 2 = 56.7 p.s/f 
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were held constant while eg varied. Figure 14 shows that T 
increases slightly from its minimum value of 263.6 revolutions as 
c B 3 i ncreases - This plot as well as the curve Stabekis obtained are 
shown, in Fig. 14, whereas the out-of-plane graph is given in Fig. 

15. 

The rotational damping constants were next varied, resulting in 
the graphs of Figs. 16 and 17. Again as the end body damping 
increases, the time constant of the least damped mode decreases. The 
values 

c Bi = 50^0 f.p/rad 

k 1 = 1000 p/f 
k 2 = 56.7 p-s/f 
were held constant. 

The time constant of the least damped mode as a function of k^', 
the nondimensional cable restoring constant, appears in Fig. 18. 

»r 

When k-j approaches zero, the equilibrium cable length approaches ■ 
infinity according to equation (32) and the definition of k]. This 
also represents an unstable situation from consideration of the nec- 
essary condition on the sign of a fi in Eq. (38). As Fig. 18 indicates, 
T also becomes large for small k^j. 

The last parameter incremented for this system was k£, the non- 
dimensional cable damping constant. From Fig. 19, the time constant 
of the least damped mode decreases only slightly as k 2 increases. 
Stabekis found that a minimum value of T occurred for k 2 approximately 
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57 p* s/f . 

The transient responses for this example with 
C£. = 5000 f «p/ rad 
kBi = 15,500 f.p-s/rad 
ki = 1000 p/f 
k 2 * 56.7 p* s/f 

are shown in Figs. 20 to 27. A comparison of the transient responses 
of the identical end mass case with the transient responses of the 
general case reveals the effect of unbalancing the space station- 
counterweight system. 

In the identical system, the amplitudes and frequencies of a-j 
and c *2 are more similar than in the unidentical case. In the uniden- 
tical case, the larger end body exhibits larger amplitude and lower 
frequency motion (Fig. 22) than the smaller end body (Fig. 23) and 
the bodies of the symmetric system (Figs. 10 and 11). 
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VI. C. Numerical Results of Gravi ty-Gradient Effects 

Equations (61) and (62) incorporating first order gravity-gradient 
effects were programmed for computer simulation. The steady-state re- 
sponse (zero initial conditions) was examined for various cases. The 
possibility of resonance in the in-plane motion was also considered. In 
all computer runs, 

R = 3.30557 x 10^f (2000 nautical miles altitude) 

G-M e = 1.407528 x 10 16 f 3 /s 2 

and Jb] - j b 2 = 1000 s1 ‘ f2 = J Ci ' I c 2 - 

Figures 28, 29, and 30 show the steady-state motion of the in-plane 
variables for the identical end mass system described in Section VI A. 

There is no out-of-plane motion for the case of zero initial conditions. 

The amplitudes of the X, and motions are in the order of 1 CT? degrees 

and for this very small amplitude motion, the response of a-j is equal to 
the response of c^. This can be seen from an examination of equations 
(61c), ( 61 d ) and (60) which verifies that for the identical system and zero 
initial conditions, the a-j and a£ responses would be expected to remain 
in phase. In considering the transient response of this system for small 
initial perturbations (Figs. 8 to 11), it can be seen that the effect of 
gravitational torques here would be negligible. 

The remainder of this section deals with the possibility of in-plane 
resonant motion due to gravi ty- gradient torques. To investigate resonance, 
the system parameters were varied, and the roots of the in-plane system 
characteristic equation, (38), examined numerically for the torque-free 
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system, until a natural frequency was found that was equal to 2e n> That 
this value of natural frequency is close to in-plane resonance is evident 
from the expression for T3 (equation 54), or T3, equation (60). 

Starting from the identical end body system, I B and 1^ were simul- 
taneously increased to a value predicted from the roots of the characteris- 
tic equation to produce resonance in one of the in-plane modes (Ig^ - Ic 3 = 
256,955 sl*f 2 ). All other parameters remained unchanged. Figs. 31, 32, 
and 33 indicate that the in-plane steady-state motion remains at about the 
same order of magnitude, compared to the motion depicted in Figs. 28, 29, 
and 30, but now the amplitudes of the higher frequency motion are increased. 

Since it was difficult (by this trial and error procedure) to find 
other sets of parameters which produced resonance for this identical system, 
the case in which an in-plane natural frequency was equal to 4e n was con- 
sidered. This situation was found to occur for I D = l r = 71,340 sl-f 2 

“3 ^3 

with all other parameters the same as in the identical case. The steady- 
state results are shown in Figs. 34, 35, and 36. It can be seen that there 
is a beat in the motion of the cable (Fig. 35). This beat frequency can be 
seen to be approximately equal to the difference between 4e n and the fre- 
quency mostly associated with the natural motion of the cable for this 
lightly coupled linear system. 

In all previous cases all parameters except the moments of inertia were 
the same as in the identical system of Ref. 6, e.g., 

P 1 = p 2 = 12 f 

= 56.7 p*s/f 

^3 = ^3 = 15,500 f*p*s/rad. 
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If 


P-| = P£ = 0.1 f 
and 

^2 = ^3 = ^3 = 0 » 

the in-plane motion is only very lightly coupled and completely undamped. 

For this system, the q^, equations can be approximated by an undamped, 

uncoupled, forced harmonic oscillator. Resonance with respect to a 2e n 

forcing frequency was predicted to occur for Ig 3 = 10*000 sl*f^ and 

c D = c„ = 10,073 f* p/rad. 
b 3 c 3 

The in-plane steady-state motion is shown in Figs. 37, 38, and 39. After 

four hundred seconds, the X and cable motion (Figs. 37 and 38) are about 

ten times larger than the X and cable motion (Figs. 28 and 29) of the 

identical system. The a-j^ motion (Fig. 39) is four orders of magnitude 

larger than the a-j^ motion of the identical system (Fig. 30). The ampli- 
tudes of all the in-plane coordinates also seem to be increasing with time. 

The effect of damping on resonance for this system is shown in Figs. 
40, 41, and 42 in which the amplitudes of the steady-state a-j , 2 motion is 
much smaller than in the previous undamped case. For this case, 

k 2 = 100 p.s/f 

kg 3 = k(; 3 = 15,500 f*p*s/rad. 

With Ig = Ic = 10,000 sl.f 2 , cg 3 = cc 3 = 16,100 f. p/rad resulted in 
resonance. From this comparison, the beneficial effect of damping in this 
resonant situation can be seen. 


57 . 


In order to induce out-of-plane motion of the system with gravity- 
gradient effects, it is necessary to have nonzero initial conditions. 

With the same initial conditions used for the identical system, i.e., 
zero initial velocities and 

X = 0 oi2 = -0.1 rad 

l - l^ = 0.48 f B-, = 0.1 rad 

= 0.1 rad $2 = 0 , 

the amplitudes of transient response in all coordinates (Figs. 8-13) were 
so much greater than that of the steady-state response for both the in- 
plane and out-of-plane motion that the effects of gravity- gradient torques 
were not apparent. However, because the eg coordinate had no initial con- 
dition and its motion is less coupled to the other equations in the pre- 
sence of gravity-gradient effects, the 63 response due to gravitational 
torques for this particular case is shown in Fig. 43. It can be seen that 
the magnitude of 02 during the first 300 second response is an order of 
magnitude less than that shown in Figs. 28 and 30. 
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VII. CONCLUSIONS 

For the three dimensional analysis of a rotating cable-connected 
space station system, the in-plane linear equations separate from the 
out-of-plane linear equations for small amplitude motion. For small 
perturbations on the cable's orientation out of the original plane of 
rotation, the system will tend to rotate in a plane inclined to the 
original plane of rotation without affecting the spin rate of the 
system. From the out-of-plane general stability criteria, positive 
damping is necessary about at least one principal axis on both end 
bodies in the plane of nominal rotation. 

From the equilibrium condition and the necessary condition for 
stability indicated by the constant coefficient of the general in- 
plane characteristic equation, the cable restoring constant must be 
greater than the value of the reduced mass of the system multiplied 
by the square of the system's inertial spin rate. 

From the necessary condition for in-plane stability, rotational 
restoring capability about an axis perpendicular to the nominal spin 
plane and on at least one end body is necessary for stability in the 
coordinates selected. For the case of identical end masses, positive 
damping and restoring torques about this same axis are necessary for 
stabili ty. 

In contrast to the general in-plane criteria for stability, for 
the special case in which the cable is attached at the center of mass 
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of the end bodies, damping and restoring effects must be provided on 
both end masses about an axis normal to the plane of rotation. 

The great difference in time constants resulting from the choice 

of the type of damping and restoring torques assumed here as compared 

6 

with those obtained in a two-dimensional analysis indicates that damp- 
ing and restoring capability proportional to angles and angular rates 
measured from the cable are better than damping and restoring with 
respect to the variational angular rates and angles ap 2 anc ^ a i j 2 as 
defined herein. 

The steady-state motion due to first order gravi ty- gradient effects 
was shown to be small and its influence on the transient response neg- 
ligible under nominal nonresonant conditions. Resonance was shown to 
occur for certain choices of system parameters. For cable attachment 
lengths which are small in comparison with the cable equilibrium length, 
the linear equations were less coupled and so the effects of resonance 
could be more easily identified. Also, damping may reduce the order of 
magnitude of the steady-state motion in a resonant situation. 

Further work on this system could involve redefining damping and 
restoring constants so that they are like those of Ref. 6 and incor- 
porating them into the equations of motion. In addition, the equations 
of motion could be rederived for the case in which each end body is 
attached to the cable at a more general point, not restricted to lie 
along a principal axis of inertia. 

A future examination of gravi ty-gradi ent effects may include a 
redevelopment of the complete gravi ty-gradi ent potential for this space 
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station- counterweight system and a more complete stability analysis 
involving Floquet theory. Effects produced by unidentical end bodies 
could also be considered. 
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APPENDIX 

i 

All nine equations of motion were derived using the approxima- 
tion that sin(qp) - qp and cos(qp) ~ 1 where qp is one of the out-of- 
plane variables (p 2 , 3p e 2 » y 2 > Y] ) * If » in addition, terms of third 
degree and higher in q D are assumed negligible, the nonlinear o 2 equa- 

r 

tion can be written as 

[ 2 (e-| + n)u - 2 (e-j +a)e-| + e-jJt] 

2 # 

+Pl[-B 3 3 2 3 2 £.+ + Pj* - (3 3 + 0)e-j3-j& - (3 3 + tt )^3^2- e l a 

+3-j 0-j p-j l + 3 3 0] B 2 3-| t - ( 9 -j + - ( e -j + q) 3 3 6 2 3 1 £]sin3 3 si ne-j 

• * <•* • * * •* 

+p 1 f! 3 2 ® ] 3 2 & — ^3^1 ^2^1^ + ( 3 3 + ^^2^1^ + 3 3 -t 

.2 

-2B30]3'j£ - (33 + ^)3 2 3l^ + 333 2 32^ + 3 2 ^ 3]- 33B 2 S]& 

2 

-(6-j + ^)3 2 B 2 £ + (0*j + 3 3 3 2 3 i &- ( 0 •j+^) B 2 32B] £- ( ® 1 +^ ) ( 3 3 +£2) 3 2 3-j £ 

-3-|^3 2 ] cos3 3 C.OS0-J 

2 m t m .. 

+Pl[33®l3]t + (33 + fi)0*jt - (33 + ft)33& - $•] 3*j ^- 3362^1 ^ — 2§3&2^1 ^ 

.2 .. .22 . .2. 

-3^£ - 23 3 B-jB 2 £ - ^ 3 ^ 2 ^ - + ^-(^ 3 +^) (^+^) ^]sinB cos9_ 

0 I 

.2 .2 2 . . 

+P}[3 2 3 2 £-3 3 B 1 3 2 £ +B 2 & +320^-320^^-3^23^+ (3 + <2)3^ 
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i 

! 

i 2 

-(& 3 +Q) 62 ^- 33 ^ 3 2 ^+(O 1 +fl)3 2 0 1 ^+(B 3 +fi)(e i +Q)Jt 

2 2 

+(6] + ^)B] 32 £ + (^1 + 0)6332* + 32 ^ 2 ^ + (33 + fl) 023 . 1 * 

+ (63 + ft) f^*] cos33sin 0 ] 

• • . 2 .. • » • • 

+P 2 [-y 3 Y 2 ' r 2 z+r 3 r 2 y l z +Y 3 £ -(y3+^)y]T] ^-(^3+^)7372^']^ 

+Y101Yl£ + Y 3 G 1 Y2Y1 * - ( 01+Q) Y] Y] * - (^l+o) Y 3 Y 2 Y 1 *]sinY3sinei 

+p 2f Y 2®1 Y 2 £ " y 3®1 y 2 y 1 £ + (Y3+^)9 1 y 2 Y^ + 2 Y 2 Y-]* + -y 3 * 

2 

-2Y3Y-|Y]^ _ (Y3+^)Y2Y^ + Y3Y2Y2* + Y2* Y]- Y3Y2Y] * 

* * • • . * 2 

- (0-j+ft) Y2Y2^ + ( 9 -j +^) Y3Y2T-J ^ + (e-j+Q)Y2Y-j * 

-( e l + fi )(Y3 + q )Y 2Y-}£ - Yl* Y2]C0SY3 C0S Q 1 

. . 2 

+p 2[ Y 3 e l Y l £+ ^ Y 3 +Q ) e l 2 -( Y 3 +Q ) Y 3 x - Y l Y l i - y 3‘ y 2 y 1 ^ “ 2Y 3 Y 2 Y 2 a 

-■Yi2.-2Y3Y 1 Y 2 )l-Y3Y2«'-(S 1 +SJ)Y3Y^-(Y3+n)(e 1 +«)i] s i n Y 3 S 1 - n 6 1 

.. ^ ^ 2 2 2 

+p 2 [ y 2 y 2 <1 - y 3 y 1 y 2 (1 + y 2 j - + Y3Y1 £ - Y20iY]A - Y 3 Y 2 Y-|(t 

+ ( Y 3 +n )Y 3 « - (Y 3 + 0 ) 6 1 a. - Y 1 ® 1 Y 2 «- - Y^Y^ + t^+njYgY J- 

+ (Y 3 +n)(0i+P.)!l + ( 0 1 +jj)y 1 y 2 ^ + (6 1 +r)y 3 y 2 « + Y 2 Y^ 
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/ ^ 

I +(t 3 + fi)Y 2 Y-|^ + (y 3 + q )y-j cosy 3 sin0 1 



n4: 


j THE VARIATIONAL FORM OF THE 0 ] EQUATION IS 

[X + 2(0' n +fl)5 + 2x5 + 2(0 n +«)55 - 2(Q n +n)9 2 0 2 ] 

+p't-S n ^2+29n s 2 s l + °1- 2 (®n +!2 )®l e r 2 (V Q )®n (3 2 6 1 +5 nh s l] sir ' 6 3Sin°i 

+Pl[20 n B 2 $ 2 - 2 ®n B 2 s 1 + 2(O n +f5 )®n 8 2 s l + 2 ^2®1 + a 1 " 2 ®nh 6 l 
-2(e n +f2)S 2 e 2 - (e n +n) 2 g 2 e 1 ]cos 3 3 cose 1 

+pJ[o 2 6 2 - 2{6 n +n)i 1 - a* . 2(e n+ n)y - 'y 1 - 2^2®! - ® 2 

■2® n h B 2 - ®n B 2 - (V fi >®n s l ' (®n +C! ) 2 - (V n ) 2|S J sin8 3 cos9 l 
+P][e 2 ^z - 2 ®n 8 l 6 2 + ®2 + ®n 6 1 “ 2 ®n 8 2 s l + 2 (8 n +njoi 1 +a 2 +2(s n +n)cj 1 5 

-e n B 2 + 2(o n +a)B23l + ( 0 n +r2 ) 2 + ( 0 n + ^' 6 + 2 ( 0 n +fi ) (3 l e 2 
. 2 

+(O n + ^) 0 n 0 2] cos$3 sin 0 ] 

' . . .2 .. . . . . . 

+P 2 [-e n ^2 Y 2‘" 20 n T 2 Y l +a 2 - 2 ( e n + ^) Y l Y r 2 ( e n +r2 ) e n Y 2 Y l +0 n Y l Y l] sinY 3 sin °I 

I 2 

+P 2 [2e n Y 2 Y 2 - 26 n Y 2 Yi + 2(® n +n )V2 Y l + 2 y 2 y 1 + °2 - 2V] Yl 


- 2(8n + S2 )y 2 y 2 - (s n + fi ) Y 2 Y il cosYj cos9-| 

j * 2 2 2 

+p 2^n Y r 2(8 n +Q)a 2 - a 2 - 2(e n +ft)a 2 6-Y-}Y] -2e n Y 2 Y] “Y]-20 n Y]Y 2 
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*®n Y 2 ' (®n + (®n + ^ * ( ®n + sinY 3 cos6 1 

+P 2[ Y 2V 2 ^S Y 2 + ^2 +5 n Y r 2 ® n :! '2 Y l +2( ®n +SJ )V“2 +2( ®n +n )“2' 5 

‘®n v 2 + 2 (® n +fi)Y 2 Y 1 + (®n +JJ)2 + (V 8 ^ 4 + 2 (V n)Y l Y 2 


* ? 

r\ . > 


+ ( e n + ^) e n Y 2^ cos ‘ Y 3 s1 ‘ n3 i = 0 


For p -j = P 2 = the linear 9^ (x) equation becomes from equation 
(28a) 


X + 2{0 n + fl)<3 = 0 


or. 


X + 2{ e n +■ ft) 5 = c = constant (A) 

where c is determined from the initial conditions. The linear l 
(5) equation is, (after employing the equilibrium condition). 


•• ^ ^ ^ ^ l i > 

5 - 2(8 n + ft)x - (0 n + Si) 6+ k-j6 + k z 6 = 0 (B) 

Solving Equation (A) for x and substituting the resulting expression 
into (B) gives 

& + k 2 5 + [k-j + 3(® n -f sij ]<$ = 2c(0 n + ft) (C) 

Now at steady-state, 6 = 5=0 so that from equation (C), the steady- 

state solution for 6 can be written 

2c(e n + 0 ) 
k] + 3(e p + ft) 2 


6 


(D) 
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Substitution of equation (D) into (A) shows that for steady- state, 
X is constant. 



4(e n + n) 2 
k] + 3(e n + si) 2 


constant (E) 


Thus, for this case, the system acquires a new spin rate and a new 
equilibrium cable length for c f 0. 
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The following two programs were used to obtain the transient 
responses: 


// XEC LIST - " ' 

// ' F OKI RAN " - - 

♦ EXTfNDEC PRECISION . 

♦ IOCSUAC3 PRINTER ,CISK ) __ _ __ 

♦ONE WORD INTEGERS 

♦LIST SOURCE PROGRAM 
REAL KI ,K2, LO 
~ ' 0 1' MEN SI OO R C 0 V M ? 6 C) 

D1 KEN'S l CM SC AL ( 18} 

_ COK/'CN MPRN r .NEXT 

COK v CM "i US), nu ( 1 6 ) , R 1 , r<> , R 1 s * R2 S , R 1 2 ;'r 1 V , R 2 V , A 1 S M , R 2 S v t R 1 2 N , 0 T H , 

1 W t lr* N , n M S t T ini N f :»J,wS»TDTH,DlHS,Xic:i,XI52fXIB3tXICl,XIC£.tXIw3 
CCK^C-’i X“Ul , <K82 » XK33 , XKC X , XKC2.XKC3, XC31, *C32 , XCB3 , XCC1 , <CC2 , 

~TxEC 37xK 1 i C £ C 0 ~ ‘ " " ' •" 

CO. i V C.N XJI , ;<J2 ,XJ3. XJ4 

__ C0KVC1N P A '•< (•' ( 5 ) ’ 

COKf'ON r. L ' ' ‘ ' ' 

EQUIVALENCE INCOf'Nt 1 ) .NPP.NT > 

DEFINE FILE 221 ( 1 , 2bC , U, I OKK Y } 

DATA SC A L / 0 . i- , 0 « t 7 . 0 , 12.0, 12.0,0. 27 ,'6 .A 0.48, 1.5, 0 . 7 , 0 .* 4 , 6 . 6 , 6 . 0 
♦,12.0,6.0,1.0,6.0,6.0/ 

PAR V U}=C.C _ 

parm'i 2 ) - 3CC: . n * ■' . 

P A R ‘i f 3 ) ~ 0 . 2 5 
PAR?-' ( 4 ) = 1 . CE-C5 

NPR NT = "" 

NEXT =1 • 

SUP =0.0 _ _ _ 

'“ tiu T 1 = 1,18 “ 

U ( I ) =0 . G 

suk=suk+ i.o/scal m_ 

rbu(Ti'=r.o/scAL('n ■ ■ ■* 

00 2 1=1,18 

_ 2 cum=oum/$i>’ _ 

' * ~ Vi = C . 050* 3 ", i 4 1 5 9 26 5 3 5 / 180. 0 

DTH= 32.0*3. 1415926535/180.0 
KN*OTH+W 
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' ■XK2*600Vb 

XK0 = XKl*X?'2/(XM+X?’2) 

l f ( fC l — y.ruv a\s ) 1 i i J 3 * .1 2 

13 ' 'WRITE (\PrfM . l*j) 

1!> FORMAT { *01 HE C A ‘3 L £ IS NOT STRONG ENOUGH FOR THE GIVEN 

_ 1 1 / * 1 1 > 

" STOP 

12 f L = ( K I*LO+X;'U* ( RHCH RHC2 ) **.NS >/ ( Kl-XKU*KNS) 


SPIN RATE 


ivRITUNPRNT , 17 )U 

17 FOR v A H AOX , * IHt tGU I L I 3 "< I U. v LENGTH IS *^1.2, /'l') 
Xl = XKU«E!.«f:L 

X IRl=piCCG.O/XI 

XI B2 = tlOOCO . C/X I 
X 1 fc 3 = & & A C G . C /X I 

X IC 1 = X 1 til 

XIC2=X1P2 

X I C 3 = X 1 0 3 


XKfi 1 -1B5CC • C/X I 
XKfi2*'i550d .C/X I 
XK l* 3 = 15500.0 /XI 

XKC1-XK01 

XKC2=XX32~ 

XKC3-XK63 


XCB l = 5 JOC . C / :< I 

’xCB2 -5 c b . 7 ;7 x i 

XC ti 3 = 5 0 00 . 0 / X I 

XCCl*XCal _ 

XCC2=XCt? 

XCC3=XC03 

XJ1 =JX! 32- X J B1 5 12. 0 
XJ? = ( X I L : 2+ X I 0 i J / 2 . C 
XJ3*(XIC2-XIC1 W2.C 


X J A = ( X ] C ? 4 X I C 1J / 2 0 

RiVRHDl/EL 


R2*KHt!2/Ejl 

R1S=R1*R1 

R2 S = R2*Jv2 
R12=R1*R2 

R1N=-R1 

R2F*-K2' 

R1S.V = -RIS 
K2S. y =-R2S 

R12M = -.R12 

XKi-Kl/xyu 
XK? = K2/X*'U 
UELG= 1 LO-EL)/fcL" 
XXK1=XKI/V.NS-1.0 
XXK2 = XK2/VmN _ 
X X K B 1 = X K ti 1 / i\ \ 
XXKB2»XKR2/hN 


X X K R 3 = X K ■) 3 / i. N 
X X K C 1 = X R C l / i% N 
XXKC2 = XKC2 / v.M 
X X K C 3 = X K C 3 / V. \ 

X X 0 3 1= X C 5 I / v;S 
XXCP2-XCf2/HNS 
XXC 3 3 = XCf'»3 / ir.NS 
XXCCUXCC1 / ivNS 
XXCC2 = XCC2/V.NS 
XXCC3 = XCC3 / V.NS 


S 
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WRI 1 EtNPRNT 
WRITE (N'PRNT 
WRl T C < \PP.M 
K I TtMRPRNT 
WRI TG ( MPRK T 
WRI T F { \ P R N T 
WRi TG (NPu-a T 
W R 1 T t ( \ P R \ T 

wr l re (:-;prm 
Wri ietf.pxM 

WRI TEUPRNT 
W’R I T t ( \PR\T 
VrI TECil'R.'iY 
V.RI IE (<\PK\T 
WRl TS (\PRNT 

writfInprm 

WRITE- (NPRNT 


?4 ) HI 
21 ) 

25) KlUl 

26 ) X l ft 2 
2 7 ) X I ft 3 

21 ) 

26 > XKIJl 

2 9 ) X n i 2 

3 G ) X K t> J 
21 ) 

'31 ) XCdI 
3 2 ) X C c ? 
3 2 ) X C ft 3 


, x [Cl 
, X I C 2 
, X 1 C 3 

, XKC1 

, xKC? 

, XKC 3 

, X CC 1 ' 
» XCC2 
f XCC3 


21 ) 

3 '. ) XR1 , XK 2 



wri rc imprnt 

21 1 





WRITE (NPRR'T 

35 1 XX Kit l 

, XXKC1 




WRI TE ( \PP,f:T 

36) X X K L* 2 , X X K C 2 




WRITE INP TXT 

37 ) X X K ;3 3 

t XXK.C3 




W R I 1 C l -■* r' R 'A T 

21 ) 





WRIT £ { ft ft R X T 

38 >XXCiJl 

,XXCC1 




WRITE! 2RNT 

39) X <C3 2 

, XXCC2 




WRI EEI 'jPRM 

< x C ri 3 

t XXCC3 




hRI T~l>J?ft\T 

21 ) 





K P. I T r 1 N P R N T 

4 1 ) XX K) , 

X X K 2 




WRI TEINPRKT 

Z?) (01 IT 

, 1=1, IS) 



21 

FORMAT ( / 1 





22 

FORMAT ( 4 fUC, 

THE PUT 

TAG CON ft 1 T I OX'S ARE 1 

Z/9X.9F 12. 

23 

FOR 2 AT ( 52 X . 

W - * 

t F 1 A * 7 > 1 • ) A f 

•UTH = 

■ , F 1 2 . 7 ) 

24 

F OR MAT { 32X, 

Rl = » 

, F 1 2 . 7 , 1 9 X » 

» R2 = 

■ ,f 12.7) 

25 

FORMAT ( 32 X, 

XI Bl - « 

i F 1 2 . 7 , 1 9 X , 

■XICl = 

• ,F1 2. 7 ) 

26 

FORMAT ( 32X, 

x i e 7 = » 

, F 1 2 . 7 » 1 9 X » 

•>:iC2 = 

‘ , F 1 ? , 7 } 

27 

FORMAT (32X, 

X I ft 3 = * 

t F-l2. 7 , 1 9 Z , 

• X I C 3 = 

' ,F12.7 1 

28 

FORMAT { 5 2 X » 

X K B 1 - * 

, F12. 7, 1 9 X , 

• XKC1 - 

• .F12.7) 

29 

FOR PAT ( 32 X, 

X K ft 2 • 

, F 1 2 . 7 , 1 5 X , 

* X K C 2 = 

•,112.71 

30 

FORMA 1 { 32X, 

XK ft 3 = * 

,F12.7,19X, 

•XKC 3 * 

• , F 1 2* 7 ) 

31 

FORMA U 32 X t 

XCBl = * 

,F 12. 7, 19X, 

•XCC1 = 

',1-12.7) 

32 

FORMAT { 32 X , 

XC ft 2 = * 

, FI 2. 7, 19X, 

• XCC2 = 

• , F 1 2 . 1 ) 

33 

FORMAT ( 3ZX, 

XC ft 3 = * 

,E12. 5, 15X, 

• XCC3 = 

♦ , e. 1 2 . 5 ) 

34 

FOR v AT { 3 2 X » 

XK 1 = * 

, F 1 2 . 7 , 1 9 X , 

1 X K 2 

• ,F12.7) 

35 

FORMA I rUX t 

XX K i> 1 = 

1 ,F12. 7, 1?X 

, > X XKC 1 

= • , F 1 2 . 7 ) 

36 

FORMAT ( 31 X , 

X X K ft 2 = 

' ,F12.7, IPX 

, »XXKC2 

= * , F 1 2 • 7 ) 

37 

FORMAT 1 31 X , 

XX K ft 3 = 

* , F ) 2 . 7 , 1 ft X 

, » X.XKC2 

= * , F 1 2 . 7 ) 

38 

F 0 R v A I ( 3 i X t 

XXL ft L = 

• ,1 12.7,1 £X 

, ‘ XX CC 1 

= ' , F 1 2 . 7 J 

39 

FORMAT (3 IX, 

XXCft 2 = 

• t F 1 2 . 7 f 1 e X 

, • X XCC2 

= ' , F) 2. 7 ) 

4 0 

FORMAT ( 31 X , 

X X C ft 3 = 

• ,E12.5, lex 

, * XXCC3 

- ' , E12. 5) 

' “ ““ ~4T 

FOR w AT( 3 1 X , 

XXKl 

* 1 F 1 2 . 7 , 1 6 X 

, 1 XXK2 

= • , F 1 2 . 7 ) 


WRI TGI 221 * 1 JMCCftN 





CALI EXIT 






EMC 






XEC 

L c 5 ( 22 


L 1 

1 ,SSCCP) 



UJ 4 
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// jns 


OOOF COO f * COOb 


OOOA 


LOG URlVli 
CCCC 
CCC t 

CUC2 

V2 f 11 


CART SPEC 

occr 
00 c A 

CCCb 


cart avail 
ccof 

_ C 09 A 
COOd 


PHY DRIVE 
0001 

_ 0002 ... 
'0000 


ACTUAL 161, CONFIG 1 CK 


// xe; list 

// CLP 

♦OfLtrh’"" "'SSK2A * ■" 

// FORTRAN 

*eXTb'iOb'C precisi CM 

* I C!C St 1 A C 3 PRI\I ER ,‘C ISK j "" 

*ONE kOKU INTEGERS 
_jLJ ST SOURCE PRCCkAK 

F X T E fiNAL $ S !. *■! * is'SG L T 
D I KENS I ON hCRK ( “ , 16 1 , NCOKM ( 260) 

CO K MON _ N ? R \ T , \ t X T 

COKMCN U( 13 ) • DU ( L6 ) » R 1 , R2 « R 1 S7 R2 S VrT?\ Y[ v . R2 V , ?.l E> t F 
1 * 1 1 k V , N S , T h t T c. t Vi S » T 0 T h t C T « S » < t B 1 t X IR 2 r X I £ 3 i R I C 1 , X I C 

COKKCN XK3 l , XK£2 , XR23 , XKC1 ,XKC_2 , XKC3, XCB l , XC22 , <Ct*3 , 

1 XCC 3 iX ; 1 * < K 2 , C E L C, ' 

CCJKVC'i XJ1, XJ2, \J3,XJA 

COWC'J 2ARM(-j> 

COP . V C \ £ L 
COKKGN T,C(9,9] 

COKKCN YY ; 36 J 

SOU I VALENCE (NCC-"’\{ i ) , NR ANT ) " 

ECU 1 VAL EN’Cfi ( C { l , 1 ) , C 11 ) , I C ( l , 2 ) , C 1 2 ) , < C ( 1 , 3 ) , C 1 3 > , { C 

1 .CIS) , tC( 1 ,6 1 ,C 1 fc) , IC( 1 , 7) ,CJ 7) , |C( t ,8) ,Clf 1 , (Cl 

2, 1 ),C2J ), (C 12,21 ,C22 1 , ( C < 2 3) ,C23) , 1 C { 2 , A ) t C2 A ) , ( C ( 7: 
3 6 I , C 2 6 ) » ( C ( 2 , 7 ) , C 2 7 ) , < C f l , c ) , C ?. <\ ) , ( C ( 2 , 9 ) , C 2 9 > 

COUI VALENCE t C < 3 , 1 > , C31 ) , (C (3,2) ,C32) , IC l 3, 3) , C33 ) , ( 

1 l 3, S ) , C 33 ) , ( C ( 3 , 6 ) , C 36 ! . ( C ( 3 , 7 ) , C3 7 > , i C ( 3, 6 ) , C33 ) , { C 

2A. 1) ,C4l ) , <Ul4,2> ,042) , IC{4 , 3) ,C43) , ( C ( 4 , ) ,C44) , iC( 

3 ' & j t C4 6 ) , ( C ( A , 7 ) , C A 7 ) , ( C ( A , & ) t C A B ) , ( C ( A , 9 ) , C4 9 ) 

tall VALchCI- <C { Y lY, C01J , ; C ( 'V, 2 ) ,CS 2 ) , (C{ 3, 3 ) ,C53 ) , ( 

1 C 5 ,6 ) , C56 1 , < C { 5 . / ) , C 5 7 ) , { C ( i> , e > , CO 2 ) , < C ( 5 , 9 > , CS9 ) f ( C 

2 6 1 2 ) t C 6 2 1 , ( C E 6 , 3 ) , C 6 3 ) , ( C ( 6 , 4 ) , C 6 4 ) , ( C ( 6 , 5 ) , C 6 0 J f ( C ( 

3 , 7 ) » C 6 7 ) , ( C l : . o ) , C t - ) , < C ( b , 9 ) , C 69 ) , ( C ( 5 , 5 ) , C 5 0 } 
EQUIVALENCE ( C ( 7 , 1) , C 7 I ! , < C ( 7 , 2 ) , C 7 2 ) , ( C ( 7 , 3 ) , C 7 3 ) , ( 
1 * 7 ' 5 > . C 7 5 > » (C [ 7, 6) , C76) , (C ( 7 , 7 ) ,C77 J , ( C( 7 , B) ,C78 ) » (C 

2 E , 1 ) , C 3 1 ) , ( C i R , 2 ) , U ; 2 ) , ( C ( E; , 3 ) t C i; 3 ) , ( C ( f> , A ) , C bA ) , ( C ( 

3,6) , CU6 ) , I C ( a , 7 ) ,CS7 ) , ( C ( ii , o ) , CE3 ) . ( C < P . , 9 ) , C89 ) 

FCtl VALENCE ) ,CV1 > , IC (0,2) ,C92) , (C (9, 3) , C93) , (' 

1 ( 9 i 6 ) , C 9 S ) , (C(9,6),C96), { C ( a , 7 ) , C 9 7 ) , ( C ( 9 , 8 ) , C9 6 ) , ( C 
DEFINE FILE 22 1< l , 2 60 , U , I D. v : ; . Y ) , 220 ( 1 2<M , 30 , U , NEXT ) 

_ F.E-ADl 22 1 ' 1 JNCCP'J 

LRITE( >, 12) NEXT 

12 FORK AT ( IX, I j ) 

_ CO 2 1 = 1,9 

UO' .2 ' J» 1' “,9 " ' 

2 C ( I , J)=G.C 
C. 1 1 = 1 . C 


S V , 912 V . 0 TH , 
, X I C 3 
CC 1 , XCC ^ , 


1 , A ) , C 1 6 ) , (C { 
, 9 ) , C 1 9 ) , ( C ( 2 
51, C 20), ( C ( 2 , 

(3, A ) , C 3 A > , (C 

3.9) ,C 39) , (C( 
» ii ) ,CA5 J , (C (4 

(5.4) ,'COA ) , ( C 
6, 1 ) ,C61 ) , (C( 

, 6 ) , C 6 6 ) , < C ( 6 

(7.4 ) ,C74 ) , (C 

7.9) ,C79 ) , i C < 

, 0) ,CF>5) , (C(& 

(9.4) ,C94), (C 

9.9) ,C99). 



C 1 3 = R1 
C1A=R2 

__ C 22 = l . 0 .. 

C 3 1 = K 1 

C33-H1S+XIH3 

C 3A = R 1 2 

~ C 4 1 = V J , 2 ‘ 

C 4 3 = K 1 2 

_C4A=R?S+XIC3 

C55-1 , C 

C66=XJ2 

C 77- X J2 __ 

COd = X JA 
C99- X JA 

_ Vv R l T C { N PP.N T » A 0 ) 

ao" r (iR’y a r ( * l * ) 

CALL RKC’SX { PARK , IJ , CU . 13. I HLF ♦ SSCjNP » SSOUT tWORKJ 

_V: R I TE: ( » 5C )_I HlF 

~ ' 50 " F 0 RKAT ( • 0" I H IF = f t 13 > * ‘ ~ 

ir;R l TE ( 22 l 1 1 ) WCCf-'.N 

CALL FXIT 

' c K r 
// DL'P 

❖ STORLCI *S UA SSX2A 1 OOOB 

* F I L £ S (22 C V S f>£0 w ) » ( 2 2 0 , 3 $ F I V ) 



